Interactions of cosmological gravitational waves and magnetic fields 
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The energy momentum tensor of a magnetic field always contains a spin-2 component in its 
anisotropic stress and therefore generates gravitational waves. It has been argued in the literature 
(Caprini & Durrer [Jl) that this gravitational wave production can be very strong and that back- 
reaction cannot be neglected. On the other hand, a gravitational wave background does affect the 
evolution of magnetic fields. It has also been argued (Tsagas et al. 0|, [S:]) that this can lead to a 
very strong amplification of a primordial magnetic field. In this paper we revisit these claims and 
study back reaction to second order. 
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I. INTRODUCTION 



Wherever we can measure them in the Universe, mag- 
netic fields of 0.5 to several micro Gauss are present. 
They have been found in ordinary galaxies Q hke ours, 
but also in galaxies at relatively high redshift Q and in 
galaxy clusters [11 . It is still unknown where these cosmo- 
logical magnetic fields come from. Are they primordial, 
i.e. generated in the early universe 0], or did they form 
later on by some non-linear aspect of structure formation, 
like the Harrison mechanism which works once vorticity 
or, more generically, turbulence has developed [^? 

In addition, once initial fields are generated, it is still 
unclear whether they are strongly amplified by a dynamo 
mechanism or only moderately by contraction. Since the 
cosmic plasma is an excellent conductor, the magneto- 
hydrodynamic (MHD) approximation can be employed 
which implies that the magnetic field lines are frozen 
in during structure formation. Therefore, as long as 
non-linear magnetic field generation can be neglected, 
the magnetic field scales inversely proportional to the 
area, so that B/p^/^ is roughly constant during struc- 
ture formation. Here p is the energy (or matter) density 
of the cosmic plasma. For galaxies, with a density of 
about pgai ^ lO^p this means that simple contraction 
will enhance magnetic fields by approximately 10'^, p is 
the mean density. Hence, if no dynamo is active dur- 
ing galaxy formation, initial fields of i?in ~ lO^^Gauss 
are needed. On the other hand, non-linear dynamo ac- 
tion can enhance the magnetic field exponentially by 
a factor up to 10^*^, so that initial fields as tiny as 
Bin ^ 10^^^ Gauss might suffice However, since this 
enhancement is exponentially sensitive to the age of the 
Universe, it remains unclear how it can generate the mag- 
netic fields in galaxies at redshifts of 2: ~ 1 or more, where 
the age of the Universe was at most half its present value 



reducing the amplification factor to less than 10^. 

Another problem of cosmic magnetic fields is that pri- 
mordial generation of fields usually leads to a very blue 
magnetic field energy spectrum. 
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where M — 2 for "causally" produced magnetic fields [l3| 
and M ~ for typical inflationary production mecha- 
nisms |11| . Such blue magnetic field spectra are strongly 
constrained by their gravity wave production [l[ and can- 
not lead to the large scale fields observed today. The only 
solution to the problem might lie in an " inverse cascade" 
by which energy is transferred from small to larger scales. 
Since within the linearized approximation each Fourier 
mode evolves independently, such a cascade is inherently 
non-linear. Within standard MHD is has been shown 
that only helical magnetic fields can lead to inverse cas- 
cade [9]. 

In this work, we want to address a weakly non-linear 
effect which has not been considered in 0, namely the 
interaction of gravitational waves and magnetic fields. 
We shall study how this interaction can modify the mag- 
netic field spectrum. We also re-interpret a finding by 
Tsagas [^], where the interaction between gravitational 
waves and magnetic fields has been interpreted as "res- 
onant amplification" . Similar conclusions are drawn in 
Refs. [13, However, in this last article it is noted 

that the amplification can take place only on super- 
horizon scales. And even though Ref. [l3| does men- 
tion that there is no amplification in the long-wavelength 
limit, they do not really quantify this statement. 

We show that the build up of magnetic fields due to 
their interaction with gravitational waves is at most log- 
arithmic and thus comparable to the generation of grav- 
itational waves by magnetic fields. 

Furthermore, in Ref. [l^ it has also been pointed out 
that the super- horizon "amplification" is independent of 
the fact whether the plasma is highly conducting or not. 
This seems physically reasonable as currents generated 
by electromagnetic fields can act only causally, i.e. on 
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sub-horizon scales. An animated discussion on this sub- 
ject foUwed the above pubheations and can be found in 
Refs. [l^, Here the role of a finite conductivity in 
an expanding Universe is addressed but controversal final 
conclusions have been reached. 

The main advantage of our treatment is that we ex- 
press the relevant results entirely in terms of physical, 
measurable quantities, which renders the interpretation 
straight forward. We actually find for the density param- 
eters of second order perturbations that, once the scales 
considered are inside the horizon, 
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II. THE BASIC EQUATIONS 

We work in the MHD approximation, where we as- 
sume high conductivity. The electric field is then small 
compared to the magnetic field in the baryon rest-frame 
which we take to be the frame of our "fundamental ob- 
server". In addition, we assume the velocity u^^ of this 
fundamental observer to be vorticity-free and we ne- 
glect acceleration. According to Frobenius' theorem u is 
hyper-surface orthogonal and we can choose spatial coor- 
dinates in the three-space orthogonal to u. Furthermore, 
in the early Universe which is of interest to us, the dom- 
inant radiation and baryons are tightly coupled so that 
the energy flux is also given by u and we can set the heat 
flux q ^ 0. In our vorticity free frame, the magnetic part 
of the Weyl tensor, Hij is related to the shear simply by 



curler, 



as one probably would expect naively. Even though most 
parts of this result can already be found in the above cited 
papers, they are interpreted there in a different way, and 
especially in Eq. ([2]) it is not always noted that the factor 
^GW s-lways has to remain small. 

The paper is organized as follows. In the next section 
we set up the fully non-linear equations for the evolution 
of magnetic flelds in the relativistic MHD approximation. 
We use the 3+1 formalism and closely follow the deriva- 
tion given in Ref. [l6j . Since we are mainly interested 
in gravitational waves, we specialize to the vorticity free 
case. In Section HI we consider linear perturbations. We 
solve the linear perturbation equations for gravitational 
waves and magnetic fields with given initial conditions. 
We also derive the evolution of the corresponding energy 
densities. This part is not new and mainly included for 
completeness and to fix the notation for the subsequent 
Section IV, where we solve the second order equations. 
On this level the gravitational waves interact with the 
magnetic field. We calculate the second order magnetic 
field generated by this interaction and show that for rea- 
sonable values for the first order perturbations, its energy 
density remains always much smaller than the energy 
density of the first order contributions. In this sense, 
one cannot speak of resonant amplification. In Section V 
we summarize our results and draw some conclusions. 

Throughout this work we use the metric signature 
(— , +, Conformal time is denoted by t and we ne- 

glect the background curvature of the Universe, K = 0. 
Spacetime indices are denoted by lower case Greek let- 
ters, fijiy, while lower case Latin letters, i,j are used for 
spatial indices. Most of our calculations are performed 
in the radiation dominated era and we shall often use the 
expression 

a{t) = H-.^afj, 



where curl is the 3-dim curl on the hyper-surface normal 
to u. Here a is the shear of u given by 



e 



and 



Note that the normalization of u implies — u'^u,^.^ cx 
uq:^. The gravito- magnetic interaction can then be de- 
scribed by the following equations, see 
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-D ai-j - K-V„Hjj - gQ^Hy 

-3cr(,;"i?j)„ - -KCr(i"Hj)„ , 
1 „ 2 

-ie^-iJp+3p+i^B^)-2.^ 



(4) 
(5) 
(6) 
(7) 



Here Eij is the electric part of the Weyl tensor, p and p 
are the energy density and pressure of the cosmic fluid 
which is assumed to follow the motion of the baryons 
(like, e.g. radiation before decoupling), k = SttG is the 
gravitational coupling constant and Bi is the magnetic 
fleld. We have neglected the electric field in the above 
equations, since we assumed it to be much smaller that 
the magnetic field, i.e. B^ ^ . The covariant deriva- 
tive in direction u is denoted by V„ and the brackets 
indicate symmetrization and trace subtraction. 
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1 



for the scale factor. 



is the Laplace operator on the hyper-surface orthogo- 
nal to u. The scalars B^ and are simply ct^ = aijcr^^ /2 
and B^ = BiB\ 
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In addition to this we have the Einstein equation, the 
spatial part of which yields 



7^^ 



+ ^kI1ij - -^ecTy + cr„(jCr"j) , (8) 



and its trace, the generalized Friedmann equation, 

le^ + ln = Kp+^KB^ + a^ . (9) 
o Z on 

Here TZij is the Ricci tensor on the spatial hyper-surface 
and TZ is its trace. 

From this system we can derive second order equations 
for an and Bi which are 
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find that the behaviour in time of the induced second 

(2) 

order magnetic field Bl is the same on super-horizon 
scales (up to uncertain logarithmic corrections). We in- 
terpret this as the insensitivity of super-horizon pertur- 
bations to plasma properties like conductivity. 

Inside the horizon, we neglect the source term. This 
is motivated by the Green function of the damped wave 
equation obtained when linearizing (jlip . which rapidly 
oscillates on sub-horizon scales. For Eq. ^ it is not 
the Green function but the source term al]^Bj,^ which 
oscillates when fci ^ 1, since gravity waves start oscillat- 
ing at horizon crossing. Therefore again, the sub-horizon 
amplification is unimportant. The same conclusion is ac- 
tually drawn in Ref. ^13i] . where the fluid velocities are 
not neglected. 

In the following we shall consider these equations in 
first and second perturbative orders with respect to a 
spatially flat Friedmann background, 



We neglect a possible spatial curvature of the background 
and work with conformal time t. The time dependence of 
the scale factor a is determined by the Friedmann equa- 
tion, 



-pa 



and 



p — —3(1 + w)p 



and 



V„V„B, - D^B 
'1 

-up- np 
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B, = 



<J^,VuB' + 2ea,,B^ + 2{Vu<y^,)B' 
3 

--nH^jB^ + a"^(Tj)nB^ + curlJ^ . 



(11) 



Eq. ([TT|) can be obtained from Eq. (40) of [T^ when set- 
ting Ai = 0, LUij — and qi — 0. Ji stands for the 
3-dimensional current. Eq. (jlip is obtained without ne- 
glecting the electric field. The term curlE^^, which is 
present in the original Maxwell equation which reduces to 
Eq. ^ ii Ei = Q 16] results in the Laplacian term D^Bi 
and terms proportional to the wavenumber k times the 
electric field [see Eq. (40) of [l^]. We have neglected 
these latter contributions in the above equation, since 
they are only relevant inside the horizon {kt 3> 1), where 
we can neglect the source term of the equations, as we 
shall argue in the following. 

In a regime of low conductivity we can neglect also the 
current in Eq. pT|) and the magnetic field obeys to the 
above wave equation, while in a very hight conductivity 
case we should directly set the electric field Ei = from 
the beginning and solve Eq. ([6]), obtaining a power-low 
behaviour with respect to time for Bi. In both cases we 



III. FIRST ORDER PERTURBATIONS 

A. Magnetic fields 

A background Friedmann Universe can of course not 
contain a magnetic field since the latter always generates 
anisotropic stresses Ily 7^ which break isotropy. When 
considering a magnetic field as a first order perturbation, 
Eq. ^ leads in first order to 



(12) 



For this we use that to lowest order u — a^^dt and 
{yuB)i — a^^{dt — a/a)Bi. Furthermore pij — gij — 
a?5ij and 8 = ia/a?. This is solved by 

B^ (x,t) - B^M—, 



5^(1) (x,t) ^ If^i^ 



(13) 



The average energy density of the first order magnetic 
field is then given by 



1 



(in) ^ ^/^4(i) 



(14) 
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Here, we assume that the first order magnetic field has 
been generated by some random process. Hence B^lj^ is a 
random variable and (• • •) denotes the expectation value. 
We assume also that this random process is statistically 
homogeneous so that (p^'') is independent of position. 

B. Gravitational waves 

For the gravity wave equation we consider a Fourier 
component 

(7y^'=(x,t) = a«(k,i)Q„(k)exp(ik-x) , 



with solution 



Here Qij(k.) is a transverse traceless polarization ten- 
sor. We assume that the gravity waves are statisti- 
cally isotropic and parity invariant so that both po- 
larizations have the same averaged square amplitudes. 
For the amplitude cr*^^-'(k, t) we obtain to first order the 
usual tensor perturbation propagation equation (neglect- 
ing anisotropic stresses of the cosmic fluid) 



= 



(15) 



where Ji = a/a denotes the co-moving Hubble param- 
eter, Ti — aH, where H is the physical Hubble pa- 
rameter. We now rewrite this equation in terms of the 
dimensionless variable S(i)(k, t) = cr(i) (k, t)/(a?j6) = 
(T(i)(k, t)/(3-ffa?j). We have normalized by the factor 
1/a?^ in order for the quantity S to be independent of 
the normalization of the scale factor. This is not true 
for a which is cr (x af^. In this way, S can be directly 
related to observable quantities which are of course inde- 
pendent of the normalization of the scale factor. Equiva- 
lently, we will make use of the variable B that is defined 
as i3 = ^/KB/{3Hain) in order to be independent of the 
normalization of the scale factor, as well as E. In terms 
of S the above equation becomes 



S(i) -3(H-w)HS(i 



(1) 
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= 0. (16) 



In the matter or radiation era, the solutions to this lin- 
ear homogeneous differential equation are well known in 
terms of Bessel functions. We are mainly interested in 
the radiation epoch, where w = 1/3. During radiation 
domination the Universe expands like a{t) cx t such that 
H — Ha = 1/t. We can therefore express the scale factor 
as 

a{t) = Hina^t . (17) 
In the radiation dominated Universe Eq. (fTB]) reduces to 

S(i)-^S(i)+(fc2 + |)s(i)=0, (18) 



cx {ktf [JAkt) + yiikt)] 



(19) 



where jn and j/„ denote the spherical Bessel functions of 
index n [ll| . 

We distinguish the super- and sub-horizon behaviors. 
In the long wavelengths limit, z = fct ^ 1, we have 

^4 



lini2:^ji(z) ~ — 

2— >0 6 



lim z yi{z) cx —z 

z— tO 



Taking into account only the faster growing mode, we 
obtain 



or equivalently 

S(i)(*) - S(i) 



kt 



fci < 1 



kt < 1 



(20) 



(21) 



The quantity directly related to gravity waves is however 
given by cr(i) = 3-ffa?jS(i) , for which we obtain on super- 
horizon scales 



^(i)(0 - ^(1) — 



kt-^l 



(22) 



A direct consequence of this is that the "gravity wave en- 
ergy density" is constant in time outside the horizon, as 
we show in the next sub-section. Of course the notion of 
"gravity wave energy density" and "gravity wave" is not 
strictly well defined for wavelengths larger than the size 
of the Hubble horizon. We shall just use the expression 
which is valid inside the horizon and call this the "grav- 
ity wave energy density" by analogy. It has a physical 
interpretation as a true energy density only once it enters 
the horizon. However, whenever this quantity becomes 
of the order of the background energy density, we know 
that perturbations become large and we can no longer 
trust linear perturbation theory. 

Let us also consider the short wavelengths limit where 
fcf 3> 1. In this limit we can approximate 



^(i){t)^{kt) 



, cos(fct) 
cos(l) 



S(i)(fct= 1) 



/ct > 1 , (23) 



where the initial constant S(i)(fct = 1) stands for the 
value of S(i) when it enters the horizon and can be ob- 
tained from Eq. 



Y.(^^){kt = 1) ~ S 



The behavior of gravity waves on sub-horizon scales, 
kt':S> 1, is then given by 



a(i)(i) ~ 3alH{ktf^^^E^,){kt 



1) 



(24) 



We shall see that in this case the gravity waves energy 
density decreases like 1/a^, as it has to be for true gravity 
waves which are massless modes. 
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C. Energy Densities 

As a first physically important quantity, let us discuss 
the energy densities of these first order perturbations and 
the corresponding density parameters. 

The magnetic energy density is 



with Eq. (fT4)l . this becomes 



(1) 



Stt 



(25) 



For the magnetic field density parameter at scale k this 
yields 



dlogfc ~ 3(27r) 



dlogfc 



(30) 



Let us now consider gravity waves. The gravity wave 
energy density in real space is given by 



P^GW 



8ttG a2 



(31) 



(26) 



In the radiation dominated universe under consideration, 
the density parameter of the first order magnetic field is 
therefore given by 



(1) 



Pc 



The density parameter f2^^ is constant in time. Both, 
the background radiation and the magnetic field which 
is frozen in, scale in the same way with the expansion 
of the Universe. As long as the magnetic field density 
parameter fi^^ is much smaller than 1, the magnetic field 
can be considered a small perturbation. 

This is the result for a constant magnetic field. We 
also want to consider a stochastic magnetic field. In this 
case B(a;) is a random variable and its spectrum is given 

by Q 



gb^ 



(l)in 



3 m 



= n 



(1) 



a2(i)B(x,t) 



(27r)3 



d3fcB(k)e 



(B.(k)i3;(q)) - i27TrS{k-<i)V.,{k)vi'l{k) . (29) 

Here the basic time evolution of the magnetic field cx a^^ 
has been removed so that, to first order B(k) is indepen- 
dent of time. VijCk) — Sij — k~'^kikj is the projection 
tensor onto the plane normal to k. The tensorial form of 
the spectrum is dictated by statistical isotropy which also 
requires that Vg]^ depends only on the absolute value 
k = |k|, and by the fact that B is divergence free. The 
Dirac delta is a consequence of statistical homogeneity-'^ . 
In this case we obtain 



1 



(27r)687r 
1 f dfc 
"fc 



d^fc / d3Q(B(k)B(q))e'^ (''-') 



(27r)3 



^k^-Pi'lik) 



dfc dp^^)(fc) 
fc d log fc 



where the factor l/a^ comes from the fact that the dot 
denotes the derivative with respect to conformal time 
and the difference of a factor 4 in the normalization as 
compared e.g. to [l7| comes from our definition of the 
perturbation variable [gij — a'^{6ij + 2ft,y )]. In Eq. (pij) 



hij is considered as tensor field with respect to the spa- 
tial metric da so that there are no scale factors involved 



(27) in raising or lowering indices, hi 



hi'' — h'-K For sim- 
plicity we shall keep this convention is this section for all 
spatial tensors. 

To lowest order the shear is given by cr|j^ = ah ij. Fur- 
thermore, the fact that cr^j'' is transverse and traceless 
together with statistical isotropy determines entirely the 
tensor structure of the power spectrum. 



(q)> = 



(27r)^5(k-q)X,,,„(k)7'(Vn(fc) 



where [l| 



(^^) Mijlmik) = SilSjni + SimS-jl - SijSim + k '^{Sijkik,-, 



+ 



in. J 



Silkjkfji Sim^ikj Sjikik-f^ 



Sjmkiki) -\- k kikjkikjji 



(32) 



We have M^^ij — 4, which takes into account the two 
polarization degrees of freedom. Therefore, considering 
that also for the shear we do not multiply by the scale 
factor while raising or lowering indices, (Xy = cr*-', we can 
write the gravity waves energy density in terms of cr.y as 



Pgw 



8nG a4 



(33) 



For the contribution to the energy density per logarith- 
mic frequency interval we then obtain 



dlogfc 



(2^)3G 

18 
(2^)3G 



k^ri'Hk,t) 



k^V^^\k,t) 



1 



( — 

a 



(34) 



^ One could also add a term which is odd under parity but wc 
disregard this possibility in this work |20ll . 



where we have used the relation Ui 



alently 'Pi^\k,t) = 9H^afj^^''{k,t). Finally, we can 



2^4 -d(1)/ 



SHaf^Yiij or equiv- 
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write the gravity wave density parameter as 

dr!g^(fc,i) _ 1 dpg]^ 



dlogfc 



Pc d log k 



(27r)3 



(35) 



We have now to distinguish between super- and sub- 
horizon modes. Using our super-horizon resuh for = 
cr(i)/(3ifa?„) where fci < 1 



we find 
dlogfc 



18 



(•in 



(27r)3G 

18 
(27r)3G 



tt2 '^PgW in 



dlogfc 



(36) 



For the last equal sign we made use of Eq. (fTT]) . Hence on 
super-horizon scales the "gravity wave energy density" 
is time independent. Then, of course the gravity wave 
density parameter grows like a^, 



di7g^(fc,t) drigl;" 



d log fc 



dlogfc 



fci < 1 



where 



dlogfc 



487r 
(2^ 



(37) 



(38) 



Inside the horizon, fct ^ 1, we have to insert the expres- 
sion of S(i) given by Eq. ((23)) in Eq. pi]) , which yields 



dlog/c 



(2^)3G I 
1 



Hi 



[kUnY 



/Oin 

V a 



(39) 



For the density parameter we obtain in a radiation dom- 
inated background 



dn'^aUk) ^ 247r 
dlogfc 



1 



(27r)3 

1 / 1 Van 

2 \kt^J "dlogfc ' 



(l)i„ 



fct > 1 . (40) 



Inside the horizon, the gravity wave density parameter 
is constant in time as is natural in a radiation domi- 
nated Universe. Note that this agrees, up to the factor 
1/2 which comes from averaging cos^(fci), with Eq. fHT]) 
at horizon entry, where (a/ajn)* = (fctin)"^. Large scale 
gravity waves from inflation, are "amplified" for a long 



time before entering the horizon, i.e. they have ktm ^ 1. 



Only if 



df}g|^(fc) /dlogfc 



is small for all values of fc, 



perturbation theory is justified. Therefore it is not suf- 



ficient if 
that (fctin) 



di7g^7diogfc 



(l)in, 



is small, but we actually need 
be small. This is better 



GW /dlogfc 

understood if we write the energy density in terms of the 
metric perturbation. In a radiation dominated Universe 
the " growing" (not decaying) mode solution for the met- 
ric perturbation is 

hij{k,t) = eij(k)hinja{kt) , 

where (k) is transverse traceless and Jq is the spherical 
Bessel function of order 0. Using Jq ~ — ji and Eq. ([3T|) 
yields 

(1) ^,^3 k'{\hn.m{kt) 

With Pc - 3i?V(87rG) = S/iSTrGa^f) and (|/ii„|2) 
we find 

dflGW 

dlogfc 



3[ikty]t{kt)]k-'Ph ~ iiktft'Vh , if fct < 1 



(41) 

Hence if the metric perturbations are small for all values 
of fc, i.e. k^Vh ^ 1 this impHes 

dflGW 



dlogfc 



« (kty 



df}g|,^7diogfc 



< 1 



Therefore the requirement (fctin) 
is equivalent to the requirement that the metric perturba- 
tions be small on super horizon scales [note that jo{z) ~ 1 
for z < 1]. 

Before we go to the second order, let us stress this 
point once more, because it is the origin of the confusion 
in the literature. Inflation generates gravitational waves 
with an amplitude 

k^Vh^ f ^]'<io-i" 



where Alp is the Planck mass and iJinf denotes the scale 
factor during inflation. The maximum value of lO"^'^ is 
the maximum tensor fluctuation from inflation allowed 
by the cosmic microwave background anisotropics. 

However, the density parameter on super-horizon scale 
is given by, see Eq. (f4T|) 



dn 



(1) 

GW 



d log fc 



(ktY 



fci < 1 



This equation is correct for any power law background, 
a oc f^, also for matter and even for inflation. Only at 
horizon crossing, can the density parameter become of 
the order 10^^°. Inside the horizon it stays constant 
if the background is radiation. Hence Eq. (1^0]) can be 
written as 



d»gl(fc) 
dlogfc 



inf 



Mp 



fct > 1 



(42) 
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IV. SECOND ORDER PERTURBATIONS 

In this section we include all terms of second order 
in the perturbations, and we shall insert our first order 
results for them; i.e. in terms of the form (JijB^ we in- 
sert ajj'' B^j^-^ or for Ily we insert the first order mag- 
netic fields, n[]^ = ^B,^^^ - (l/3)p|-^B(i)2 ^ Eqs. H 



and pT|) . We obtain the following differential equa- 
tions for the evolution of the second order perturbations 
Bf^(x,i) andag'(x,i): 



iT^ '\/ +2B(T^ 'B^ +2V (7^ 'B^ 



.(1), 



+ 



(i:)2)(i)B(i) +curlJ, 



(43) 



(2) 



(2) 



+ ie2(l-3u;)a(f 



(1) 



l2^(l)Jl) 



.(1) 



(44) 



Taking into account that V„-Bi(i 



gether with V„ng'' : 
be simplified to 

VuVuBP - d^bP 

"4 



(1) 
t(i) 



-(2/3)ei?, 



(1) 



to- 



(4/3)en^V Eqs. 1121), (glD can 



-9V B^^-' 
3^ v„ijj 



^(1) 



(I?2)(i)b« 



(45) 



(2) 



ie^(l-3^)ag) 



(1) 



0O'(i(l)'"cr 



(1) 



2^(1)^(1) 



+2a 



,(1) 



(i?-)-'(T, 
(1) 



(46) 



We have also neglected the term curlJ^ in Eq. (|43p . Since 
it is proportional to k in the Fourier space, its contribu- 
tion is important only on sub-horizon scales, where we 
anyway neglect the source part. Outside the horizon, 
fcf <C 1, it is negligible. 



A. The second order magnetic field from gravity 
waves and a constant magnetic field 

For simplicity, and to gain intuition, we first consider 
a constant first order magnetic field. 



i3«(x,t) = B-^, 



(1) Oin 



B^,'\k,t) = Blti^6-\k). 

In this case, the convolution of -B(i) and (T(i) into which 
the products in ordinary space transform under Fourier 



(l)«in c3/ 



transformation, become normal products and the second 

(2) 

order magnetic field B^ has the same wavelength as the 
first order gravity wave which generates it. 

-4. 



Remembering that 
obtains 



oc a 



V. 



(1) 



a a 



4^(1) 



1 



2^(1)5(1)^ 



(47) 



In principle, one has to consider the corrections to 
the orthogonal spatially projected covariant derivative 
(£)2-j(i)^(i) ^YiQ tensor perturbations hij in the 

metric tensor g^^^,. Computing these corrections, they 
turn out to be equal to zero, since the magnetic field is 
transverse. This remains valid even if B^-*^^ is not con- 
stant. 

Considering the expansion-normalized dimensionless 

f2) f2) 

variable Br' = V^Bf V(eain), we obtain 



^ - nil + 3w)B, 



(2) 



fc2 



.(1) 



B 



1) /Oil 

V a 



(48) 



We investigate the behavior of the second order per- 
turbation in the radiation dominated phase. 

Moreover, since the source /^(k, t) and therefore also 

(2) 

Bl (k, i) are random variables, we want to determine 
their spectra. The first order gravity wave spectrum is 



(E«'"(k)E,y^'"(q)) = {2nrM.,MS-\k-ci)V^^ik) 



(1) 



(E«'"(k)S«j^„(q)) = 4(2^)^5^(k- q)7'^^4(fc) , 

where Mijim is the gravity waves polarization tensor de- 
fined in Eq. ((32)) . It can also be expressed in terms of the 
projection tensor Vij{k), Mijim = VuVjm +VimVji - 
'PijVim- Actually {l/2)Mi/™ is the projection tensor 
onto the two transverse traceless modes of a rank 2 sym- 
metric tensor. The power spectrum of the second order 
magnetic field i3(2) is of the form 

{BP{k,t)Bf\p,t)) = i2nfV.,{k)S'{k-p)r^^\k,t) . 

(49) 



5(1) 



(2) 

We obtain the solution for B^ (k, t) with the help of 
Green function method. 



(k,o = £ dt'git,t',k)f,{k,t' 



)• 



(50) 



Here Q is the Green function of the second order lin- 

(2) 

ear differential operator acting on B^ which depends on 
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the cosmological background. It can be determined in 
terms of the homogeneous solutions which in the radia- 
tion dominated era are simply spherical Bessel functions 
and powers. More precisely, in terms oi z = kt, Eq. 
in the radiation dominated case, w = 1/3, becomes 



)" - f-B^>' + + ^ j - fc"V.(^,k) , (51) 

where the prime denotes a derivative w.r.t. z. Two ho- 
mogenous solutions to this equation are Pi{z) = z'^jQ{z) 
and P2{z) = z-^yo{z). Defining the Wronskian, W{z) — 
P[{z)P2{z) — Pi{z)P2{z) = z^, a possible Green function 
is 

gfe.M.)^ ^-"'>^'M-;;-'"^-"'' . (52, 

The solution obtained by integrating with this Green 

(2) 

function satisfies the initial condition Bl (zi„,k) = 

(2,)/ 

B\ ' (zin, k) = 0. Any other solution can be obtained by 
adding a homogeneous solution to this one. We discuss 
the physically correct choice of initial conditions in more 
detail in the Appendix |^ For the magnetic field, the 
initial conditions chosen with this Green function seem 
adequate to us. We can now write the magnetic field 
spectrum as 



j(2)' 



j(2) 



g(z,z',k)g*(z,z",p)(/,(k,z')/;(P,^")> 



(53) 



We solve Eq. ([5T|) . distinguishing the sub- and super- 
horizon regimes. In the long wavelength limit, kt = z 
1, we have to insert the solution obtained for gravity 
waves E-j^ on super-horizon scales and given in Eq. ((20)l . 
Therefore, the source term fi (k, t) reads 



(k)B(i)'" (i?i„ain)'t' ,(54) 



and equivalently for /j*(q, t"). The power spectrum of fi 
can then be written as 

(/,(k,z')/;(p,^")> = 

16«Pnk)P](q)(Si'j'"(k)Sri'''"(k)) 

= i2TTfS\k - p)V^J{■k)h{z', z", k) . (55) 

For the function h{z',z",k) we obtain the following ex- 
pression 

h{z',z",k)^F{k)g{z')g{z") , 
F{k)^nBlyy^l{k)k-\ 
g{z') = AHlflir^z' . 



The solution for the power spectrum of the second order 
perturbation of the magnetic field can then be written as 



{Bf\k,t)B]'^\p,t)) - (2^)^P,,(k)J^(k-p) 

n 2 



dz'g(z,z',k)^/F(fc)g(z') 



(56) 



The square [• ■ ■ ]^ is simple the power spectrum V^\k, t) 
which we want to determine. Of course, the integrals 
in the square bracket are solutions to our magnetic field 



Eq. ([5T|) with source a/ F{k)g{z). Hence \j'P'^^ satisfies 
the equation 



2 

P" - -P' 



\p\^^vf{Kt), 



P=-z 

fc3 ' 



= kt 



(57) 



a^iHla^^^^fiBlyy^'lik) . 

Solving the above equation with the Wronskian method 
in the regime z = fct <C 1, one finds 



z = /ct < 1 . 



This yields 

k'V^^\k,t) 



4k 



Din 2 

^(1) 



kt < 1 



(58) 



This is the second order magnetic field power spectrum 
induced by the presence of a first order field and a grav- 
itational wave. It is the growth cx t^ of this induced field 
which has been interpreted in Refs. 0, |^ [l2| as strong 
amplification. But before drawing such conclusions, we 
want to compare the energy density parameter of JS^^^ 
with the one of a^^^ and B^^^ inside the horizon, where 
these quantities have a simple physical interpretation. 

Inside the horizon, kt ^ 1, we can no longer use the 
above simple approximation for the source term. The 
solution of Eq. (j57p with a generic source term, 



P^^\k,z) 



Vj^\k,z) 



V^i\k,z) 



S{k,z), (59) 



can be written as 

\jv'i\k,z) = j dz'S{k,z')g{z,z\k). (60) 



But, once the gravity waves enter the horizon, the source 
and the Green function start oscillating and the con- 
tribution to the above integral becomes negligible. We 
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therefore neglect the source inside the horizon and sim- 
ply match the solution at horizon crossing with the ho- 
mogeneous solutions of Eq. (|57ll given above, that are 
Pi{z) = z'^jQ{z) and P2{z) = z'^yo(z) [z = kt). Consid- 
ering the limit z ^ 1, this yields 



2k 



Din 2 

H? 

Ill 

1 

(fctin)4 



fci > 1 



The energy density 



(61) 



Cosmic Microwavwe Background (CMB) anisotropics 
on large scales (integrated Sachs- Wolfe effect). Simi- 
lar arguments yield flg^ < 10^^ on large scales (see, 
e.g. plj, H^l). Therefore, even though we agree with 
the calculation in Ref. 0], we do not agree with the in- 
terpretation. If the gravitational wave energy density 
is as small as required by the measurements of CMB 
anisotropics, flg^ always remains smaller than 17^' . Fur- 
thermore, up to logarithmic corrections, -B*-^^ inherits the 
spectrum of the first order gravity waves. 

In the next section we show that this conclusion per- 
sists also if we allow for a stochastic magnetic field. Just 
the computation becomes more involved. 



To analyze this amplification which happens mainly 
on super-horizon scales, let us compare energy densities 
after horizon entry. The energy density of our second 
order magnetic field is 



dlogfc 



1 



(2^)3 

1 

(2^ 



(fc,t) 



1 



PT'f (fc,t) 



^(^^)\(62) 
K \ a / 



The factor 1/a^ comes from the fact that we have to 

in order to compute the 
energy density, while is due to the definition of 



raise one index of (B^'s?') 



(2) 



(1) 

Bl '/air, that we gave above. The density parameter for 
B^"^^ then reads 



dlog k 



(2^) 



k^vf{k,t) 



(63) 



With H = Hinaf^/a^ we find that even though 'P^^\k,t) 
is growing like on super-horizon scales, the density pa- 
rameter grows like ^^q^- After horizon entry, this growth 

(2) . 

stops and fl^ remains constant. Inserting the solutions 

(2), 



and (Ell) for k^V^^'{k,t) gives 



An''g\k,t) 
dlogfc 



6- 



Gw(fc'^)o(l) 



dlogfc 



(64) 



on super- and sub-horizon scales. 

Hence, even though the second order magnetic field 
B(2) is growing considerably, this reflects only the growth 

of the unphysical density parameter J^qIv on super- 
horizon scales. Once this is factored in, the mag- 
netic field density parameter is not. The values for 



both. 



An 



(l)in^ 
GW 



(fc)/diogfc (fcii„)-4= dr!g^(fc)/diogfc 



and Vl'^g are at most of the order of 10 ^ and 
smaller. For the gravity waves, we have seen that 
dr2Q^"(fc) /dlogfc (fciin)^'' is just the square ampli- 
tude of the metric perturbations on super horizon scales, 
which has to be k^Ph ^ 10^^° in order not to overproduce 



B. The second order magnetic field from gravity 
waves and a stocliastic magnetic field 

In the case in which the first order magnetic field is not 
spatially constant, all the products S^-j^ (x, i)i?^^^ (x, i) 
become convolutions in Fourier space 



d3xe*'^-Eg)(x,i)5^,)(x,i) 



^^'."(k) j d3gE«(q,i)i?^-^^(k-q,t), 

where the projector Vi" = (Sf — fc^fc" projects onto the 
transverse modes. The result of this convolution is a 
magnetic field and therefore transverse. Hence this pro- 
jector is not strictly necessary. But as we shall see, it 
simplifies the calculations. 

Our equations are written in terms of the dimensionless 
expansion- normalized variables (x, t) and (x, t) , 
and we want to express their power spectra in terms 
of the power spectra of the first order random vari- 
ables (x, t) and S.^j-* (x, t) for which we assume simple 
power laws. 



i?f)(k,t) = Bj^)(k) 



a 

.3 



i3(\)(k,t) = S""(k)^ 



af„(B«'"(k)i3;^^^'"(q)) = (2^)37^,, (k)53(k - q)7^),\/Jfc) 



(1) 

in 



(1) 



^'"(k)i3r^,„(q)) = 2{2nr5\k-ci)V\,'!jk) , 



5(1) 



,(1) 



(fc) 



[B2^^;„A3](Afc)^^ forfc<fcrf, 
for fc > kd, 



(65) 



where kd is the damping scale which we assume to be 
always much smaller than the Hubble scale. The scale A 
is arbitrary, e.g., the scale at which we want to calculate 
the magnetic field. With this normalization B"^^ is sim- 
ply the amplitude of the magnetic field at scale A at time 
tin . At any other moment, the magnetic field at scale A 
is given by B'^^afja^{t). 
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Equivalently we have for the gravity wave power spec- 
trum 



(E(f'^(k)S;W'"(q)) = 

(2^)3X,,,„(k)<53(k-q)P^^4(fc), 



(1) 



(66) 



Here the transfer function T{k,t) keeps track of the de- 
terministic time-dependence of the gravity waves. In the 
previous section we have derived the well known behav- 
ior of the gravity wave transfer function which oscillates 
on sub-horizon scales, fci 3> 1, and behaves like a power 
law on super-horizon scales. For the radiation dominated 
case, 



r(fc,t) 



kt < 1 



(67) 



Starting from Eq. , we can write the following evo- 
lution equation for the second order perturbation 

(x, t) + 2nBf^ (x, t) - a^D^B^^^ (x, t) + 



(68) 



Replacing i?!^^ = 3HainB, 
we obtain 



j(2) 



K and = 3i?a?„sg\ 



2v^aina-B^i-,(x,t) 



(1)/ 



(69) 



This is the same differential equation as for the constant 
magnetic field. In Fourier space this equation becomes 



^ (k, f) - (1 + 3w)nBl'> (k, t) + BY> (k, t) X 



(2), 



/.(k,t), (70) 



where the source fi (k, t) is now given by a convolution 
2 



/.(k,t) 



(27r)3 



/Kaina7'i'"(k) x 



j d3qS«(q, i)i?^^j(k - q, t) - ^(1 + u;)7^x 
Jd^q4'^{ci,t)Bl^{k--,l,t) . (71) 

In terms of the variable z = kt we obtain again 
Eq. (|5T|) . As in the previous section we solve it with the 



Green function method. Therefore, the power spectrum 

(2) 

of Bl is given by 



(Bl'\k,t)B;''>ip,t)) = (2^)^5^(k-p) X 
{^ij^kikj)T'^\k,t) , 



-(2) 



with 



r^^\k,t) = / dz' dz"g{z,z',k) X 

g*(z,z",p)(/,(k,z')/;(p,^")> , 

where z = kt. In the radiation dominated epoch {w 
1/3) the source term reads 

/,(k,t') = ^^a,,Mt')'P,'-{k) X 

1 d3gi],«(q,i')i?ri)(k-q.^')- 
2n{t') J d3gEW(q,i')Sa)(i^-q'*') 



(27r)3 ^ a2(i') 
1 d'qi:ill"\ci)fiq,t')Bi^y^-ik-cD- 

2n{t') I d39EW/"(q)T(g, t')s4!^ '"(k-q) 



(72) 

and equivalently for fj{p,t"). To determine the power 
spectrum of fi we assume that the magnetic field -B(i) 
and gravity waves (T(i) are uncorrelated, so that 



(/.(k,o/;(p,i")> = 



16k 



H{t')H{t'- 



a{t')a{t" 



{2TTf 

d3Wd3.(I](r(q)E:W'"(s))x 



n'^(k)7',"(p) X 



(BW'"(k-q)i3rW'"(p-s)) 

= (2^)353(k_p)7',,(k)/l(t',t",fc) . 

The function h{t',t", k) is given by [l| 



(73) 



h{t',t",k) 



8k 



n{t')H{t") 



a{t')a{t") 



I(k), 



(27r)3 

I{k)^ j A\{l+^^){l+a^)v'il{q)V'iU\^-ci\), 

(74) 

where a = k ■ (k ^ q) and 7 = fc • We neglect the 
angular dependence of (1 + 7^) and (1 + a^) and simply 
set 

(H-72)(H-a2) ~ 1 . 
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We then have to solve the following integral, 

f 1/ inax(t',i") 



dqq 



A+2 





M/2 



Here we evaluate the integral only up to the scale q which 
enters the horizon at the later of the two times. All 
scales q < l/max(t', i") are super- horizon from tin to 
max(i',i"). A soon as q enters the horizon, the gravity 
wave transfer function begins to oscillate and the contri- 
bution to the integral becomes negligible. The integral 
over ^ can be evaluated; for M ^ — 2 it yields 



2 + M (1)'" (1)'"' 

l/max(t',i")j^ 

kq 



A+M+6 



{\k + q 



M+2 



\k-q\ 



M+2\ 



We shall not treat the case M = — 2, where the angular 
integral introduces a logarithmic dependence on q, sepa- 
rately. This corresponds to approximating log(A;/g) ~ 1. 
We approximate these integrals by their dominant con- 
tribution. 

• If the spectra are sufficiently red such that A + 
Af -t- 3 < 0, the result is dominated by the region 
k < l/max(t',t") and we obtain 

I{k) 167r^l^,^Bl^,^\' X 

^^''^ U + 3 A + M + 2, 

• On the other hand, if the spectra are blue such that 
A-\-M+?> > 0, the integral is dominated by its value 
at the upper boundary. 



167rS^i) 
A 



,A^ 



1 



(l)in-- ^. 



M + 3 



max(t',t") 



If, as in the previous sub-section, we can write the 
function h{t',t", k) in the form 



hit',t",k)^Fik)git')gin 



(75) 



we can proceed as we did before to obtain the results 
([58| and (|6ip . A source where the time dependence of 
the unequal time correlator factorizes is called "totally 
coherent". In the totally coherent case, the power spec- 
trum is simply the square of the solution which has as 
its source the square root of the power spectrum of the 
source [23j . In most cases, the unequal time correlator is 
more complicated than this, but the totally coherent ap- 
proximation is often quite reasonable 23]. If the source 
is totally coherent, the square root of the power spectrum 

(2) 

■pg simply satisfies the same evolution equation as S(2) 
with source term ^/Fg. 



• If A + M + 3 < 0, we can write 

F{k) = i^(fcA)^+*^+3A3 X 



{2n) 
1 



1 



9{t') 



A+3 A+M 

-B(l)inS(i) 



3^ 



• For A 

9{t') - 



- Af + 3 > 0, we set 
1287rK 1 



(27r)3 A + M 

-S(l)inS(i 



y^A+M+Q 



This corresponds to replacing 



1 



max(t',i") 



(A+Af+3) 



by 



1 



(A+A/+3)/2 



{A+M+3)/2 



which is of course not entirely correct and we ex- 
pect this to over estimate the true result somewhat. 
However, within the accuracy of our approxima- 
tions this is sufficient. To obtain a more accu- 
rate result we would have to expand the function 
h{k,t',t") in eigenfunctions with respect to convo- 
lution in time, as it is done in Ref. [23|. 

Within this totally coherent approximation we can now 
solve the problem like in the previous sub-section. In the 
case ^-|-A/+3 < we find on super- horizon scales, where 
the source is active 



ik,t) 



32ttk 




(27r)3 


HI 




6 


(-) 


, fci < 1 


Vain/ 





(76) 



On sub-horizon scales, performing the matching at hori- 
zon crossing, we obtain 



ik,t) 



k'vi'lik) 



{2n) 



Hfn 
1 



{kU 



> 1 



(77) 



If ^ -I- Af -f 3 > 0, we analyze in in more detail only 
the case A ~ —3 and Af — 2. The spectral index A = 
—3 correspond to a scale invariant gravity wave power 
spectrum as it is obtained in slow- roll inflation . The 
index Af = 2 characterizes a causal magnetic field -B(i). 
In this case, we have to solve the differential equation. 



P" - -P' 



1 



P 



a 

fc2 ' 



(78) 



a=(ai„ffin)"i?L'^4l^ 



647r 
(2^ 



kA5 
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where the source is constant in time. Detailed comments 
about the initial conditions chosen for the solution of the 
above equation can be found in Appendix El Finally, we 
can write the solution for P{z) in the case where z — 
kt as 



z < 1 



The power spectrum of B^^'' on super-horizon scales is 
therefore given by 



(27r)3 



HI 



{ktii, 



■logM- 



fc^T'W (fc)^ X 
/ci < 1 . (79) 



On sub- horizon scales, z = kt^ 1, we match the super- 
horizon solution at horizon crossing with the homoge- 
neous solution of Eq. ((75)) . as we did above, obtaining 



k^v'i\k,t) 



327rK 



(27r)3 



(ktinY 



■ log^ (fcti, 



/ct > 1 . (80) 



1. Density parameter 

Using Eq. ([62]) . we find the following expressions for the 
energy density of the stochastic second order magnetic 
field. If + M + 3 < 0, we have on super-horizon scales 



dlog/c 



1 



2887r r 
(2^ 



k'V^lik) .(81) 



This results in a density parameter for _B(2) given by 
dn'§\k,t) 



dlogfc 



1 dp^^\k,rj) 
Pc dlogfc 
df7«„(fc) df7g^(fc,i) 



6 

dlogfc 
/ci < 1 . 



dlogfc 



(82) 



Inside the horizon we obtain for the second order mag- 
netic field density parameter 



The 



dn^g\k,t) 

dlogfc 



gravity 



d log k d log k 
fct > 1 . 



(1) 



6 



(83) 



wave 



density parameter. 



dn^(^]y{k,t)/dlogk is given by Eqs. §7^) and (gD]) 



respectively. This corresponds, as in the previous section 
for a constant magnetic field, to the naively expected 



result, 



(2) 
B 



"gw"b ■ 



For blue spectra, A + M + 3 > 0, the second order 
magnetic field density parameter reads in the interesting 
case A ~ —3 and M = 2 on super-horizon scales 



dn'-^\k,t) 
dlog/c 



12 dn 



(1) 



[kty dlogfc dlogfc 



(1) 



12 



d log fc d log fc 



{ktf 



k=l/t 



log^ — 



Note that the value of 



fci < 1 



(84) 



df7g^(fc,t)/dlogfc 



on super- 



dfl')^\kt)/dl0gkt 



at hori- 



Hubble scales is affected by 

zon crossing, kt = 1/t which may well be larger than 

df7^^'(fc)/dlogfc but of course has also to be much 

smaller than 1. 

This expression grows only logarithmically faster than 

dn 



(1) 



Gwik,t) /d\og 



The growth stops at horizon entry 

where the second order magnetic field density parameter 
has acquired a factor log^(fcii„). Inside the horizon we 
obtain a density parameter of 



dn^g\k,t) 
dlogfc 



12 



dn^^lik) d»g>y(fc) 



dlogfc 

fct > 1 . 



dlogfc 



log^(fctii, 



(85) 



Up to the logarithmic correction, this corresponds to the 
result for red spectra above. 



2. Reheating and matter dominated epochs 

In order to make contact with Refs. 0, we now 
repeat the calculation in a matter dominated background 
{w = 0). We want to point out that the results we obtain 
are mathematically the same as the ones found in Q . The 
only difference lies in the interpretation. In the previous 
paragraph we have seen that, even though 



dlogfc 



1 



k'Vilik) 



and even though (fctin)^^ can become very large, this 
product is never larger than about 10"^". We believe 
that this point has been missed in Ref. !3|. 
If w = 0, the scale factor grows like a c>c so that 
Ti. = 2/t. As mentioned before, for the super horizon 
amplification the question whether the conductivity is 
high or low is not relevant. 

^From the first order perturbations, we obtain the 
same behaviour for the magnetic field B^^-* in terms of 
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the scale factor, therefore the density parameter is then 
given by 



dn^g' {k,t) 
dlogfc 



SttG 



3(27r)3 



Hi 



(86) 



The first order gravity waves on super-horizon scales now 
behaves as 



I] in 



'i„(k) — 



(87) 



Once the gravitational waves enter the horizon, they start 
oscillating and the energy density decays as radiation. 
Therefore in this case the relative density parameters for 
the first order gravity waves is 



dlogfc 



(27r)3 
kt < 1 



(88) 



On sub-horizon scales we obtain 
247r 



df^g^(fc,t) 



dlog/c 



(27r)3 
fei > 1 . 



(fc) 



(fc<in)« 



(89) 



Computing finally the induced second order magnetic 
field density parameter, we obtain the naively expected 
result on super-horizon scales 



dfl^g\k,t) 
dlog/c 



da 



d log k d log k 

for A + M ^ 



3 < 



{kty 



d log k 

(1) 



GW 



(fc,t) 



d log k 

for A + M 



X 

k=l/t 

3 > 



fci < 1 



(90) 



On sub-horizon scales the density parameter turns out to 
be given by 



dn'-g\k,t) 
dlogfc 



dn^j^\k,t) dr>^^4(fc,tfc) 

dlogfc 



(1) 



dn'-g\k,t) /Hi 



d log fc 

2 



'inf 



fci > 1, (91) 



d log fc \ Mp 

for both cases A + M + 3 < and A ~ -3, M = 2, up 
to logarithmic corrections. Here tk stands for the hori- 
zon crossing time, tk = 1/fc, and in the last ~ sign we 

have used that dilQ^(fc,tfc)/dlogfc ~ (iJinf/Mp)^ is 
the gravity waves density parameter at horizon crossing, 
which is smaller than 10"^". This means that the sec- 
ond order magnetic field does not grow larger the the 
first order one. Inside the horizon they decrease both 

(1) 

B ' 



like 



(2) 

fin stays always much smaller than 51 



as we have found in the case of a radiation dominated 
background. 



C. Second order gravity waves 



Starting from Eq. (|46p . we can write the evolution 
equation for in real space (x, t) as follows: 



•(2) 



a^^Vlf-:^7^^(l + ^).(f 



(1) 



2a(7 



1^(1) 



(92) 



The factor 1/a^ in the source part of the above equation 
comes from the fact that in Eq. (|46l) we had to add factors 
a^(i) in order to lower or rise indices. On the other hand, 
now we deal with purely spatial tensors such that aij — 
a^^ and also &ij = . 

Introducing again the dimensionless expansion- 

(2) 

normalized variable E^^^ , the previous equation can be 
written as 



Sg'-3(1 + 7«)HE,^ 



r(2) 



2w 



.(2) 



2 a 



--(1 
2^ 



2 '^h] ^ 



(1) 



n y(l) 
(■'(1) j)n 



As for the source is given by the first order pertur- 
bations magnetic field [n|j^] and the first order gravity 
waves and does e.g. not couple to the second order mag- 
netic field. Since we assume the first order magnetic field 
and gravity wave fluctuations to be independent, we can 
add the power spectra for the solutions of the individual 
source terms. 



k''V^^\k,t) = k^V'^>'\k,t) + k^V- 



(2)n, 



^3Tr.(2) GW 



(fc,t) 



where V^''^{k,t) is the power spectrum of the solution 



of Eq. ((93|) with source term II'^^ only and V- 
comes from the source terms containing E^^-*. 



(2) GW 



(fc,t) 



1. Magnetic field part of the source \k^V^^^ (k,t)^ 

Considering first the magnetic field part of the source, 
we have to solve the following differential equation in the 
momentum space (k, t) 



- 3(1 + w)n^ 



(2) 



.(2) 



k . (94) 
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where the source is given by 

2 a2 



(95) 



As before, we have to compute the unequal time correla- 
tor: 

(ng)(k,i')n:«(p,i")) = 

(2^)353(k - p)M„,„(k)/i(fc, t', i") , (96) 
where the anisotropic stresses are given by 

S«(k-q,t') . 

(l/2)A^ij'''(k) is the projector on the tensor modes. We 
have neglected a trace contribution to the magnetic field 
stress tensor since, once we project with Mi/'' , the trace 
vanishes. 

After some computation [l|, we find for the function 
h{k,t' , t") the following expression: 



h{k,t',t") 

m 



1 



1 



(87r)2 4(27r)3 



m 



1 2 



a{t')a{t") 



(97) 



d\{l + ^^){l + a^)V^^l{q)x 



^i'L(ik-qi: 



(98) 



where a = k ■ [k ~ q) and j = k ■ q. As before, we 
approximate (1 + 7^)(1 + a^) ~ 1. With this, we obtain 
the following expression for the expectation value of the 
source term: 



(/,,(k,t')/™(p,i")> = ^^^^^^¥^x 

^ ^in 

t(i) 



nl^(k,t')n;,«(p,i")) 



(99) 



(2) 

The expectation value of the stochastic variable S^^ can 
be written as 

(Eg)(k,t)E;W(p,t)) = 

{27TfS^{k-p)M,,rn{^)r^^\k,t) . (100) 



If {Il'^\]i,t')Ilrn'{p,t")) can be written as a product of 
a function of (fc, t) and {k, t"), this source is totally coher- 
ent and we can write the function h{k, t' , t") of Eq. (P7|) 
in the form 



— K 

9 



,a2(t')a2(t") 



h{k,t',t")=F{k)g{t')g{t" 



where we introduced the pre-factor of h since we finally 
need an expression for the unequal time correlator of the 
source, as in Eq. ([M)) . while the function h alone is only 



part of the correlator of the anisotropic stress, Eq. (|96|) . 
The square root of the power spectrum is then a solu- 
tion of the differential equation (IM)) with source term 
\/ F{k)g{t). Written as differential equation for the vari- 
able z = kt and setting w = 1/3, this becomes 



,(2)n 



{k,z) 



,(2)n 



(fc,z) 



,(2)n 



(fc,z) 



- Vm'-^ (101) 



As for the second order magnetic field, we distinguish 
between two cases. First we consider 2M -I- 3 > 0. The 
integral I is then dominated by the upper cutoff. The 
magnetic field is not oscillating and we therefore take 
damping scale kd as the upper cutoff. We neglect the 
slow time dependence of this scale. Using Eq. ([55]) for the 
magnetic field power spectrum, / can be approximated 

by 



Stt 



(Afcd) 



2A/4-3 



2Af + 3 

Hence the functions F{k), g{t') are given by 



m 

9{t') - 



1 



36(27r)4 2M + 3 
1 . 



,2M+3 



In the case 2M + 3 < 0, we obtain 

1 



Stt 



(Afc) 



2M+3 



1 



. M + 3 2M + 3 
This case is totally coherent and we can set 



m = 



git') 



36(27r)4 
1 



1 



M 



2M-f 3 



(Afc) 



2M+3 



1 . 



We now solve Eq. (|101[) for the two different source 
terms. 



• In the case 2M + 3 > 0, we can write Eq. (|10ip in 

the form 



4 

P" - -P' 



= kt 



P 



a 
fc2 



,(2)n 



{k,t) , 



Solving the above equation on super-horizon scales 
and following the considerations for the choice of 
initial conditions explained in Appendix [Aj we find 



P{z) 



2fc2 



z < 1 
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this gives the second order power spectrum 



so that 



36(27r)4(2M + 3) 



-\ 2 



4 / , X 2M+3 

a \ I kd 



k 



, kt^l . (102) 



k'vi'Uk) 



144(27r)4 

4 

— 1 , fci < 1 . 

air, 



(107) 



This is equivalent to a density parameter for S^^-* 
given by 



d»g|v"(fc,^) 
d log k 



d»L'!n(fc) 

dlogfc 

dlogfc 
fct < 1 



2M+3 



(103) 



Inside the horizon, the Green function oscillates 
and we can neglect the contribution from the 
source. The solution for the power spectrum is then 
given by 



k^r^^^''(k,t) 



36(27r)4(2M + 3) 



kH'lik) 



-\ 2 



Hi 
kt > 1 



4 / , X 2M+3 

a \ / kd^ 



(104) 



Therefore, the second order density parameter is 
given by the same expression, 



As in the first case, the density parameter is the 
same for fct < 1 and fct > 1, 



dn^^l^{k,t) 

d\og k 



d\og k 



(108) 



2. Gravity waves part of the source 

Let us finally consider the part of the source given by 
first order gravity waves. In this case, we can write the 
source as: 



/..(x,i)= [-|(l + 3z.)H2Eir4i 



As before, we ignore the traces that are present in the 
above products, once we evaluate them in the momen- 
tum space, since we project them out with {l/2)Mi/"^ 
afterwards. Remembering that S^- = E*^ , we have on 
super-horizon scales, where the transfer function is given 
by Eq. dSZl): 



dn^Gl^{k,t) 

dlogfc 



dn<^Blikd) 



dlogfc 
fci > 1 



(105) 



Up to logarithmic factors this result agrees with the 
findings of Ref . [l[ . 

• In the case 2M + 3 < we have again to solve the 
equation 



4 

P" - -P' 



4 \ a 



(106) 



Hence 



But now 



1 / 1 



1 



6(27r)2 y 2\M + 3 2M + 3 



^2M+3 



X 



M 



(k,t) 



1 



2(27r)3 \au 
M./"^ik) J d3gE«'"(q)I]W/"(k-q) 



2n 



(k,t) 



(27r)3 \a 



— ) M,/™(k)x 



(k,i) 



/ 



d^gE«*"(q)EW'"(k-q) 



These equations are strictly true only on super-horizon 
scales where E cx 1/a^. However, since inside the horizon 
E oscillates and the contribution from the source is neg- 
ligible, we can use this approximation. Setting w = 1/3 
we can finally write the source in the form 



/».(k,t) 



2(27r)3 

M./™(k) j d3pEif"(p)EW-(k-p) 



(110) 
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and the two-point correlation function of the source part 
reads: 

h{k,t',t") , 
1 



h{k,t',t") 



8(27r 



-^U{t',ni{k) 



81, 



ait') 



a{t") 



u{t',t") = ^'H'{t')n'{t") 



I{k) = MbdlniOi-) / d^P Mlnbf{p)M,nndf{^-p) + 



More details about the computation of h{k,t't") and of 
the four point correlation function of the gravity waves 
can be found in Appendix iBl 

Using the tensor calculus package "xAct" for Mathemat- 
ica [23 , we can compute the above products of the three 
projectors, 

M bdlm 

(k) \Mlnbf{p)M,nnd.f{k-p) + 
Mlndf{p)M„inbf{k - P) 

2(1 + +f3^ + o?0^ ~ 8af3"f + 

(Ill) 

where a = k ■ [k — p), (3 = p ■ [k — p) and j = k-p. Again 
we have approximated this angular dependence by a con- 
stant to simplify the calculations. This approximation is 
well justified within our accuracy. In order to write the 
function h{k, t' , t") ~ F{k)g{t')g{t"), we have to evaluate 
the integral / as before. We first consider the most in- 
teresting case of a scale invariant spectrum, A ~ — 3. Up 
to an infrared log-divergence which we neglect as usual 
(this divergence can be avoided if we choose A = —2.99 
instead of A = —3), we have 



BItt 
2(2^ 



a{t') 



2 1 



Therefore, the equation for y {k , t) in the radia- 

tion dominated era becomes 



P"{z) - -P' + ( 1 



z = kt , 



z 



P = \V- 



,(2)GW 



(fc,t) , 



a : 



SItt 



k^V^^lik) 



if„^ 2(2^)3 k^ 



The super-horizon solution, evaluated always with the 
help of the Wronskian method and keeping only the non- 
homogeneous part as explained in the AppendixEl is the 
given by 



z < 1 , 



that yields a contribution to the gravity wave power spec- 
trum given by 



^ ^3^(1) 

"(27r)3 
fci < 1 



12 



(112) 



For the density parameter on super-horizon scales this 
yields 



An'§^{k,t) 
dlogfc 



0.01 



An^G^{k,t) 

dlogfc 



fcf < 1 . 



(113) 



Considering now the sub- horizon limit, we obtain for the 
power spectrum the following expression: 



fci > 1 , 



(fciin)s \a. 
and the density parameter becomes 



(114) 



dlogfc 



0.02 



d»L'^(fc) 
dlogfc 



, fci » 1 . 



(115) 



On the other hand, when 2A + 3 > we have 



2(27r) 



3^(l)in'- 



2A-I-3 



a{t') 



(2A+3)/2 



In the radiation epoch the equation for ^V^^'^'^ {k,t) 
reads: 

z V z-^ 



a= JJll ^fc3/2p(i) 1 

Y2(27r)3 2A + 3 ^'"^ ^fcii„)6 • 

Solving the above equation in the long wavelengths limit, 
we find: 



P(z) ~ -z9/2-4 , 

V ; 2 



z < 1 , 
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where the exact pre-factor depends weakly on the value 
of A. For the power spectrum this results in 



SItt 



^ ~ 8(27r)3 

and inside the horizon this reads 
SItt 



12 



(116) 



3t,(2)GW 



16(27r)3 
1 



A:t > 1 , 



(117) 



Translating this to the density parameter as above, we 
obtain 



dr>g^^(fc,t) 

dlog/c 



0.2 



0.2 



dlogfc 



dr!g^^(fc,t) 

dlog/c 



dlogfc 

fc< < 1 , 

d»g^(fc) 
dlogfc 



(H) 



fc=i/t 



-2A-3 



{ktf 



fci > 1 . 



(118) 



V. SUMMARY AND CONCLUSIONS 

In this work we have studied the evolution of stochas- 
tic cosmic magnetic fields and gravity waves up to sec- 
ond order in the perturbations. We have especially cal- 
culated the density parameters of the generated second 
order perturbations. We start with density parameters 



which are 



dn^j^\k,t)/d\ogk and dn'^^l^{k,t)/dlogk 

related to the first order magnetic field and gravitational 
wave power spectra in Section [ill CI Since tensor pertur- 
bations grow on super-horizon scales, the gravity wave 
density parameter grows on super-Hubble scales and only 
becomes constant once the perturbations enter the hori- 
zon. For perturbation theory to be valid, we have of 
course to require that these density parameters are much 
smaller than unity. As we have seen in Section IIII C[ 



to require that dflQ^{k, t) /d log k 
also on sub-Hubble scales, is equivalent to 



is smaller than one 



1 



dlogfc (fcii„)4 



inf 



Afp 



< 1 



(119) 



Here we summarize the new results on the density pa- 
rameters for second order perturbation on sub-horizon 
scales. For magnetic fields, we obtain 



d»g^(fc) ^ dr>g^(fc)d»W(fc) 

dlogfc 



d log k d log fc 



tfc > 1 



(120) 



up to numerical constants and logarithms which are be- 
yond the accuracy of our approximation. Hence, it is 
not correct that the presence of gravity waves resonantly 
enhances a first order magnetic field. The second order 
density parameter is quite what we would naively expect 
and it is much smaller than the first order perturbations 
as long as the latter are small. Also on super-horizon 
scales, the second order magnetic field density parame- 
ter is always much smaller than the first order one, see 
Eqs.JUl) and (|84)) . 

Since the growth comes from super horizon scales, con- 
ductivity is not relevant for this result. We have shown 
that also in a matter dominated background we obtain 



dn'i\k) ^ dn\i^ik) 

d log k d log fc 



(1) 



inf 



kt=l 



dn's\k) 
d log fc 



dn'^'ik) 

Mp J dlogfc 

dn'^\k) 

d log k ' 



(121) 
(122) 



hence no significant amplification. 

Second order gravity waves are induced on the one 
hand by the anisotropic stresses of the first order mag- 
netic fields and on the other hand by the quadratic terms 
in the evolution equation for try which are, e.g., of the 
form (TimtT™ and similar expressions. In Section FlV C II 
we have shown that the second order contribution from 
anisotropic stresses on sub-Hubble scales is of the order 
of 



d»[|i"(fc,^) 
d log k 



d log kd 



d log k 



if 2Af + 3 > 



if 2M + 3 < 



(123) 

both on super- and sub-horizon scales. Note that the 
above expression is continuous at 2M -1-3 = 0, where 
both expressions scale like (fcA)'^ and are independent of 
kd- One should point out that we neglected the slow time 
dependence of the damping scale. Correctly one has to 
choose the value of the damping scale at horizon crossing, 
kd{tk) with tk = 1/fc. Depending on the magnetic field 
spectrum, the resulting gravity waves come mainly from 
the small scale magnetic field, if its spectrum is blue 2M+ 
3 > 0. In this case the gravity waves power spectrum is 
always proportional to fc^ . In our case of a simple power 
law magnetic field spectrum, this behavior is maintained 
for all fc < fcrf. If the magnetic field spectrum is red, 
2Af -I- 3 < 0, gravity waves depend on the field at scale fc 
and their spectrum is the square of the _B-field spectrum. 
In the first case, the non-linearity leads to a 'sweeping' 
of magnetic field power on small scales to gravitational 
wave power on larger scales. This can be regarded as an 
'inverse cascade' of small scale magnetic field power into 
large scale gravity waves. But in no case can the gravity 
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wave density parameter become larger than the one of 
the magnetic field, which has to be much smaller than 
one, for perturbation theory to be valid. 

A similar result was already obtained in Ref . 1] . Con- 
trary to this reference we have no logarithmic build-up of 
gravity waves. This comes from our different treatment; 
we directly calculate the shear cTy and not the tensor per- 
turbation of the metric, hij. In this way we loose the log 
term which corresponds to the homogeneous hij — con- 
stant solution on super-horizon scales to which we are 
not sensitive. However, in our more qualitative work, we 
do not want to insist on log terms which we neglect in 
this work also in other places. 

The second order gravity wave density parameter in- 
duced by first order gravity waves is given by 



an 



(2)S 
GW 



ik,t) 



dlogfc 



dlogfc 



A:<> 1 



(124) 



on sub-horizon scales. 

Adding both contributions we find 





2 , 


dng^(fc,t) 


d log kd 




d log k 



if 2M -I- 3 > , fci > 1 



dlog 



d log k 

if 2M 



da<!4(fc,t) 

d log k 

3 < , fci > 1 
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APPENDIX A: GENERAL SOLUTION OF A 
DIFFERENTIAL EQUATION WITH THE 
WRONSKIAN METHOD 

Here we discuss in detail the Wronskian method with 
which we find the solution of the differential equations 
in this paper. If we have a inhomogcncous linear second 
order equation with inhomogeneity S{z), its most general 
solution is of the form 

P{z) = ci{z)Pi{z) + C2{z)P2{z) + aiPi(z) + aaPzW , 

where Pi{z) and P2iz) are two (linearly independent) 
homogeneous solutions which we suppose to be known, 



W{z) = P1P2 — P1P2 is their Wronskian, and 
ci(z) = - 



, S(x) , , 



C2{z) 



dx 



W{x 
S{x 



P,{x) 



The particular solution given by the first two terms is 
such that Pinh(^) = ci{z)Pi{z) -f C2{z)P2{z) vanishes at 
z = Zin and also Pl^Yii^in) = 0. The general solution is 
obtained by adding a homogeneous solution, Phomiz) = 
aiPi{z) + a2P2{z) with arbitrary constants ai and a2. 
Let us first consider the example given in Eq. (|106p . 



P" - -P' 



1 



P 



a 

fc2 ' 



where a/k'^ is a constant source term. The homoge- 
neous solutions are given by Pi{z) — z^ji{z) and P2{z) — 
z^yi{z) and the Wronskian determinant reads 

W{z) = z^ . 

In the regime z <C 1 we can approximate the spherical 
Bessel functions by powers and we find the following gen- 
eral expression for P{z): 



a f z"^ z'^ 



Z^ ~ ZZ\, 



+ aiz + a2Z 



(Al) 

where we have used the fact that, when z <C 1, we can 
(125) approximate Pi{z) ~ z** and P2{z) ~ — z. Now it is im- 
portant to notice that the second and the fourth terms 
of the inhomogeneous solution (|Aip have the same func- 
tional behavior as homogeneous solutions and we can al- 
ways choose ai and 02 such that the homogeneous part 
cancels them. This is actually always true for the contri- 
butions from the lower boundary of the inhomogeneous 
solution. This may sound pedantic, but it is very im- 
portant in this specific case as the second term in (jAl|) 
dominates if it is present. In our analysis we have always 
subtracted such "homogeneous contributions" and only 
kept the "minimal part" , which in this case is 



z < 1 



(A2) 



This procedure is important and it is responsible for the 
results which we have obtained. Wc justify it also by the 
fact that the first order solution has exactly the the same 
time evolution as the homogeneous term and therefore a 
term cx z^ present at early times, should be included in 
the first order perturbations. Once the wave number has 
entered the horizon, z 3> 1, the Green function starts to 
oscillate and the additional contribution to the integral 
can be neglected. We then can match the inhomogcncous 
solution at horizon crossing to the homogenous one at 
later times. Up to matching details which we have not 
considered, this yields 



P(z) 



2k 



^z cosz , z 2> 1 . 



(A3) 
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In the same way, we deal with Eq. ([78] 



(A4) 



The homogeneous sohitions are Pi{z) — z^jo(z) ~ 
and P2iz) = z'^yo{z) ~ —z. These approximations are 
vahd for z <C 1. Using again the Wronskian method, 
we obtain the following general solution on super-Hubble 
scales, z <C 1: 



P{z) = ^ ( Z^logl — ) -Z^+ZZin) +aiz-'+a2Z , 



z < 1 , 



(A5) 



Here, the homogeneous solution parts are — z^ and 
zzin, therefore we can identify the solution due to the 
presence of the source again as 



P(z)~ —z' log l^—j , 2«1. (A6) 

On sub-horizon scales this becomes, up to matching de- 
tails which only modify the phase and have an irrelevant 
effect on the pre-factors. 



P{z) ~ — log(fciin)z COS z , z > 1 . (A7) 

If the source term depends on z, the details of the cal- 
culation as well as the results change somewhat, but the 
basic argumentation remains the same. We therefore do 
not repeat the z-dependent examples which arise in this 
work here. 



APPENDIX B: THE FOUR-POINT 
CORRELATOR OF GRAVITY WAVES 



Gaussian, therefore we can apply Wick's theorem. The 
we can write the products of four gravity waves S^^^ as 

'"(p)S«;"(k - p)S:« '"(s)E;« '"(q - s)) = 
(E« '"(p)I]:« '"(s)) (I]W;"(k - p)S;« '"(q - s)) + 
(SW '"(p)I];« '"(q - s)) (EW '"(k - p)S,t^ ''^(s)) + 
'"(p)S« ''^(k - p)) (S;}^) '"(s)E;« '"(q - s)) . 

(B3) 

Once the double integration is performed, the last term 
contributes a constant oc 6'^ (k) which can be disregarded 
(a background term). Integrating the remaining two 
terms over (Ps, we can eliminate one of the two 6- 
functions which come from the expression of the two 
point gravity wave correlator. Using the reality con- 
dition, E*j(k) = Eij(— k), and the expression for the 
two-point correlation function of gravity waves given in 
Eq. ([66]) . we then obtain 



J d'p J d3s(EW'"(p)EW'"(k-p)S:?^'"(s) X 



-Mlndf{p)Mmnbf{'k - p) 



(B4) 



The above equation is symmetric in k and q, as well as 
under the exchange of the first and second pairs of in- 
dices. Moreover, it is symmetric under the exchange of 
the first index with the second and the third with the 
fourth. This suggests us to write the two point correla- 
tion function of the source term as 



Starting from Eq. piOp , we compute the two-point cor- 
relation function of the source term (/y(k,t')/*^ (p,t")), 
which is given by 



(/,,(k,i')/;.(q,i")> 



U{t',t")M^/'^{k) X 



i2ny 



Mj'ici) I d'p I d3.(EW'"(p)sW'"(k-p) 



E:«"'(s)I]7;^"'(q-s)) 



fd 



(Bl) 



where the function U{t',t") contains the all time- 
dependence of the above expression: 



a{t') 



a{t") 



(B2) 



To compute the four-point correlator, we assume that 
the random variables that describe gravity waves are 



(/,,(k,t')/;,(q,t")) = {27TfS'{k~q)M,,rc{k)x 



h{k,t',t" 



(B5) 



since the tensor Mijrc bas the same symmetries. 

To obtain an expression for the function h{k, t', t"), it is 
sufhcient to calculate the trace of the above two point cor- 
relator. We hence should multiply the r.h.s. of the above 
equation and of Eq. (|Bip by A^'-^'"'^(k). Then, setting 
them to be equal and remembering that Ai^^^"^ Mijrc = 8 
[20| . we obtain 

8(27r)353(k - q)/i(fc, t") = Uit', i")(5^(k - q) x 



X"''™(k)X'"*„(q) j d?p \Minbf{-p)M„,ndf{k-p)+ 

Mlndf{p)M.nnM{^)\ ^Eb)^Sil(|k - Pl) 

(B6) 



20 



with 

h{k,t'X 



:U{t',t")M hdlra 

(k)x 



8(27r)3 

/ d^pV^^l{p)v'il{\k-v\)^ 

Mlnbf{p)MmndfO^ - P) + A^tod/ (p) A^mnfc/ (k - p) 

(B7) 



Finally, we have to perform the above product of three 
polarization tensors, defined as in Eq. ([5^ . To achieve 



this aim, we use the free source package "xAct" for Math- 
ematica 25|: it is sufficient to define a three dimen- 



sional flat metric and the projection tensor 'Pij(k) = 
Sij — k~'^kikj onto the plane normal to k. Then, we 
can express Mijimi^) in terms of this projector as 



M 



ijlm 



Defining the angles between the three directions as a = 
k ■ (k — p), (3 = p ■ (k — p) and 7 = fc • p, we obtain the 
expression given in Eq. (jf If p . 
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